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Abstract 

We calculate the classical gluon field of a fast projectile passing through 
a dense medium. We show that this allows us to calculate both the ini- 
tial state gluon production in proton-nucleus collisions and the final state 
gluon radiation off a hard parton produced in nucleus-nucleus collisions. 
This unified description of these two phenomena makes the relation be- 
tween the saturation scale Qs and the transport coefficient q more trans- 
parent. Also, we discuss the validity of the eikonal approximation for 
gluon propagation inside the nucleus in proton-nucleus collisions at RHIC 
energy. 



LPT-Orsay 06-41 

1 Introduction 



The advent of a new generation of colliders, RHIC and LHC, has stimulated 
the development of new tools for the understanding of high energy and high 
density systems. The purpose of studying heavy ion collisions at high energy is 
to create such dense systems, leading eventually to evidences for the existence 
of the Quark-Gluon-Plasma (QGP) predicted by the QCD phase diagram. Very 
early, it appeared that one has to distinguish between final state interactions, 
occurring after the hard parton production, and initial state interactions respon- 
sible for hard parton production in nucleus-nucleus collisions; this motivated the 
study of deuteron-gold collisions at RHIC where final state interaction are ab- 
sent. 

1 Unite mixte de recherche du CNRS (UMR 8627) 
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The physics of proton-nucleus colhsions at high energy turned out to be very 
rich in new features compared to proton-proton coUisions, for instance: high 
Pt suppression at forward rapidities, Cronin enhancement (also observed at low 
energy), centrality dependence of spectra, etc. [3,4]. These results have been 
confronted with the theory of the Color Glass Condensate (CGC) which de- 
scribes the nuclear wave function at high energy [1,2]. Basically, it extends 
small coupling QCD calculations to a region, the so-called saturation regime, 
characterized by the hard scale Qs (saturation scale), where high density effects 
do not allow one to apply the usual perturbative QCD. 

Historically the idea of saturation of the gluon distribution at high energy has 
been introduced in the early 80 's [5-7] as a necessary condition for unitarity. 
Indeed at high energy gluons dominate the dynamics and the growth of the 
number of gluons, driven by the BFKL evolution equation, violates unitarity. 
This equation resums large log(s) effects but disregards gluon recombination 
[8,9]. A first version of the CGC known as the McLerran-Venugopalan model 
pointed out the interest of the semi-classical picture for the description of high 
density systems. [10-12]. A couple of years later, a more sophisticated theory 
has been built which incorporates the main BFKL features and extends them 
to the non-linear regime: the saturation regime, leading to the BK-JIMWLK 
equation [13-24]. 

On the other hand, several works inspired by the Landau-Pomeranchuk-Migdal 
(LPM) effect in QED [38-40]: the BDMPS formahsm, based on Feynman di- 
agrams calculations [25-31], and the Zakharov formalism (Z) based on path 
integrals [32-37], have been dedicated to the study of final state interactions 
in nucleus-nucleus collisions, especially to understand the large suppression of 
observed large pt hadron spectra compared to proton-proton collisions. The 
basic idea is to explain this suppression by the energy loss of the initially pro- 
duced hard parton. It looses part of its energy by radiating gluons when passing 
through the dense medium formed after the collision. In [41-44], Wiedemann 
(W) has provided another treatment of the problem also in terms of Feynman 
diagrams. 

In the present work, we give a new and simple derivation of the radiated gluon 
spectrum, recovering the BDMPS-Z-W result. We provide a universal formu- 
lation for both gluon production in proton- nucleus collisions [45], and gluon 
radiation off a produced hard parton in nucleus-nucleus collisions. The basic 
idea is to calculate gluon radiation off a high energy projectile passing through a 
dense medium in the semi-classical picture. To do that, we solve the Yang-Mills 
equations for the radiated gluon field SA^ , which is treated as a perturbation of 
the background medium field Aq. This medium could be a nucleus (cold mat- 
ter) in proton-nucleus collisions or a hot medium produced after the collision 
of two heavy nuclei. This approach turns out to be a very useful framework, 
which avoids many technical problems, making the picture clear. Obviously, 
the obtained classical field is a function of the medium source density po (or 
equivalently, the medium background field A^). The source distribution po is 
a random quantity. Thus, to calculate observables one has to average over all 
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possible source configurations with a given statistical weight W[po]: 



(O) EE / V[pompo]Ol 



Pol- 



(1) 



The differential average number of produced gluons is given by the formula: 



where e'^.^s^ is the polarization vector of the radiated gluon, and q = (q^ = 
uj,q^ = q']_/2uj,q^) its momentum. 

In section 2, we derive the gauge field induced after the interaction of a fast 
projectile with an unspecified medium characterized by the statistical weight 
W, by solving the Yang-Mills equations. Choosing the gauge as the light cone 
gauge of the projectile allows us to give a simple derivation, and we show, in 
section 3, that it is straightforward to deduce from the classical field the induced 
radiative gluon spectrum in nucleus- nucleus collisions, and in the approximation 
of independent scattering centers we reproduce a well-known formula derived 
in the more formal approach mentioned above (see BDMPS-Z-W papers listed 
above) . 

In section 4, we consider gluon production in proton-nucleus collisions. Assum- 
ing the interaction time to be much smaller than any other time scale appearing 
in the problem (that amounts to consider gluon propagation inside the nucleus to 
be eikonal) wc recover the fci-factorization formula which rcsums high density 
effects in the nucleus [45,46,48-51]. Then wc discuss the validity of this ap- 
proximation, and show that it fails in some kinematic regime probed at RHIC 
energies. Finally, wc summarize in section 5. 

2 The gauge field of a fast moving parton (or 
iiadron) passing through a dense medium 

We consider a massless parton (or a hadron described as a collection of partons) 
moving with the velocity of light, in the direction. At some time it passes 
through a dense static medium of size L ^. The medium is described by the 
following current, in the medium rest frame. 



We are interested in the induced radiative spectrum off the hard parton inside 
the medium. Knowing the simplicity of solving the Yang-Mills equations in the 
light-cone gauge of the parton = [45], we perform a boost of velocity 
/3 ~ — 1 (the cross-section is Lorentz invariant) which affects only the medium, 

■^The assumption of a static medium (static scattering centers) in parton energy loss has 
been used in BDMPS-W-Z formalism and in the Gyulassy-Wang model [52] . 
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Po(a:^a;x)<5^°. 
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namely the parton remains in the direction and the medium is pushed very 
close to the light-cone in the x~ direction. 

In the boosted frame the medium field can be checked to be ^ [46] : 

A^ = -S-'^^po{x^,x+), (4) 

In the light-cone gauge of the fast parton, = 0, the Yang-Mills equations 
read 

- ig[A\d+A'] = J+, 
[D-,d+A-] - [D\F'-] = J-, 

d+F-' + [D-,d+A'] - [D^,F^'] = . (5) 

Wc assume that when the parton traverses the medium, it induces a perturba- 
tion 5A^ of the strong mcdimir field A^ (linear response): 

A^' = A'^ + . (6) 

Similarly for the conserved current 

= + SJ^" . (7) 

Keeping only terms in the Yang-Mills equations which arc linear in the fluctu- 
ation JA^: 

-d+{d^SA^^) = sj+, 

aSA' - 2ig[AQ,d+SA'] = ^'{^^,SA^'), 

a6A- -2ig[AQ,d+6A-] = 5J-+2ig[d'A^,SA']+d-{df,SA'') 

-tg[A^,d^6A^]. (8) 

The parton current obeys the conservation relation 

d+6J- +D-5J+ =d+5.J- +d-5J+ ~ig[AQ,dJ+] =0. (9) 

With the initial condition 5J+(x+ = to) = 6{x~)p{x±), where to is the source 
production time, namely the production time of the hard parton in nucleus- 
nucleus collisions. In writing this current, we assume that the hard projectile 
propagates in the x~^ direction of the light cone; therefore, its interaction with 
the medium is eikonal: it only gets a color precession when passing through the 
medium. The solution of (9) reads 

SJ+ = Uix+,to,x^)d{x-)p{x^)e{x+ -to) , 

SJ- = -e{x+)p{x^)S{x+ -to) , (10) 

^For any gauge choice for the medium field in its rest frame the A+ component is suppressed 
by the boost. 
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where U is a Wilson line in the adjoint representation of the gauge group: 

U {x'^ ,tQ, x±) = exp 

where 7^ denotes the time ordering of the integrals along . The current 5J~ 
corresponds to the propagation of an antiparticle moving in the opposite direc- 
tion of the hard parton we are interested in. It is necessary to get a conserved 
current. Making use of the the first equation of (8) (seen as a constraint because 
it contains no time derivative) in the last equation for 5A^ we get : 

U6A- ~2ig[A^,d+6A-] ^ 2ig[d'A^,5A']+5J--^{d-5J+-ig[A~JJ+]) , 

. 

Thanks to the current conservation (9), 6J^ cancels out in right hand side of 
equation (12). Finally, we simphfy equations (8) which reduce to 

DSA' - 2ig[A^,d+6A'] = --g^SJ+, 

as A- ~ 2ig[AQ,d+6A-] = 2ig[d'AQ ,5A'] + 2SJ- . (13) 
The first equation of (8) is solved leading to 

6A'{x) ^-J d^z0{z')G{x, z)d\U{z+,to, z^)p{z^))0{z+ - to) . (14) 

The gluon propagator G is the retarded Green's function obeying the equation 
of motion: 

(□, - 2tgiA^.T)dt)G{x, y) = 5{x - y) , (15) 
with the initial condition, the free retarded propagator 

G\x,y) = -e{x+-y+)e{x--y-)5{[x-yf)= / ^ . (16) 

The propagation of the emitted gluon, contrary to the hard projectile, is not 
necessarily eikonal. Eq. (14) has a simple diagrammatic representation, shown 
in fig. (1). The color precession of the source before the gluon emission is 
accounted for by the U that multiplies p; while the rescatterings of the gluon 
after it is emitted are hidden in the Green's function G. The component 5A~ can 
be extracted from the constraint (first equation in (8)), but it is not relevant for 
gluon production since only transverse polarizations are physical. The integral 
over z~ is restricted to the positive values: this comes from the term (l/(9+)(5J+ 
(which contains a 5{z^)) appearing in the second equation in (8) after the 
substitution of SA'^ from the constraint. 

The radiated gluon transverse field (14) is the main result of this section, it 
contains, as we will show, the physics of proton-nucleus collisions at high energy 
and of the induced radiative gluon spectrum in nucleus-nucleus collisions. 

*In light cone gauge, we can forget about the <5A^ since it plays no role in the gluon 
spectrum. 
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dz+ A^''{z+,Xi_)Ta 



(11) 
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Figure 1: A schematic representation of eq. (14). The fast projectile (thick hne) 
passes through the medium of thickness L and emits a gluon at the time z+. 
The gluon emission could also occur outside the medium, at z+ < or z+ > L. 



3 Induced radiative spectrum in A-A collisions 

We assume that the fast projectile is a hard parton produced in a nucleus- 
nucleus collision at t^, and take the origin of times when the parton enters the 
medium. In coordinate space, (for .t+ > L) the gluon field (14), amputated of 
its final free propagator, reads 



UJA\x) - -e{x-)di{JJ{x+M.xi_)p{x^)) 

- j (fze{z-)9{L - z+)e{z+ - ta)5{x+ 
X 2dtG{x,z)dl{U{z+ M\zi_)p{z^)), 



L) 



(17) 



where we used in the second term the following remarkable property of the 
Green's function (valid for j/+ < L < a;+)[46]: 



G{x,y) 



j dz~d^z±G{x,z)2d+G{z,y) 



(18) 



z+=L 



In equation (17), the first term corresponds to gluon emission occurring after 
the hard parton left the medium and the second term corresponds to gluon 
emission before the hard parton left the medium (if the parton is produced out 
side the medium, this emission could also occur before it enters the medium). 
The invariance by translation, with respect to the - variables, of this Green's 
function is due to the fact that the medium field Aq is independent of x~ , 
this can be seen in eq. (15). This invariance implies that G depends on the 
coordinates only via x^ — . 

It is useful to introduce a new Green's function, defined as 



^^(x+,a;_L;2/+,y^) = 2 J dl d+G{x+ ,x±;y+ ,yj_]l' 



y )e 



(19) 
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It is easy to verify that Gu; is a Green's function of the two-dimensional Schrodinger 
operator: 

{id' + ^+ ig{A^.T))g^{x+,xr,y+,y^) = i5{x+ - y+)6{x^ - y^) , (20) 

therefore, it can be written in terms of a path integral for a quantum particle 
of mass uj moving in a potential: 



/ 



X'r^(^) exp 



xJ7(x+,j/+;ri), 



(21) 



where r±{x'^) ~ x±_ and r"^(y+) = This path integral describes the Brown- 
ian motion of the emitted gluon in the transverse plane. Therefore, the emitted 
gluon follows a non-eikonal trajectory inside the medium ^. We recover the 



eikonal case by taking w — > cx) or equivalently by assuming x 



0: 



5^(.T+,a;^;y+,y^) = 5{xi_ - y^)9{x+ - y+)U{x+ ,y+:xi_), 
The Fourier transform of (17) gives 

d^x^e-'i^-^^ [^—U{L, to; x^)p{xj 
q + le 

[L, xr, z+,z^)dl{U{z+,to; z^)piz^))] 
The amplitude for the gluon radiation reads 



(22) 



q'SA^q) 




LO + is 






+ J d^z, 


_ [ dz+ 




J to 



(23) 



(24) 



where the gluon is taken on-shell i.e. q^ ~ (or q' = q^/2a>) . To get the 
average number of the produced gluons, we square the amplitude and sum over 
the polarization vectors with the help of the completeness identity 



(25) 



In the same way as in the CGC treatment we have to perform the average over 
the sources. First, for the fast parton, we write 



{p^{x^)p\x'^)),^ = pI{x^)5'^H{x^ - x'^) , 



(26) 



^ Since the medium is static these are the only non-cikonal corrections to gluon propagation. 
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IJ,p{x±) is the parton charge density in the transverse plane. 
From (24) (or equivalently (2)), the gluon spectrum reads 



OJ 



d{N) 





[-1 


X ( dz + 


J to 


Vuj J 


X {tvU{z+ 


z'+ 



(2^) 

.+ 

dz'+ 



(trC/(z+, 2'+, z^)dl,gl{L, y^-z'+, z'^W,g^{L, xr, z+ , z^)) 
^2fil{y^)-^ (trC/t(i, z+, y^)dig^{L, x^;z+, y^)) 

- 1) 



(27r)3q2 



j d^x±fil{x±) 



For a single parton produced at x± = 0± it reduces to 



m - 1 



<5(a5±), 



(27) 



(28) 



where i? = A for a gluon and R = F for a quark. 

In eq. (27) , x± end ?/_|_ are the coordinates of the emitted gluon in the trans- 
verse plan. The first term in (27) corresponds to the probability of producing the 
gluon inside the medium, whereas the second term corresponds to the interfer- 
ence between the amplitudes for producing the gluon outside (after the parton 
left the medium) , -first term in (23)-, and inside the medium,-second term in 
(23)-; this is illustrated in fig. 2. In the second term, note that the transverse 
coordinates of the emitted gluon and of the hard parton are the same (see the 
first term in (23) where the final gluon free propagator has been amputated). 
The last term is the probability of radiating a gluon in the vacuum, it has to 
be removed to get the medium induced gluon spectrum. This formula is quite 
general, in the sense that we do not specify the nature of the medium, so one 
still has to find a model for averaging over the medium sources. 
Now we will show that this formula leads to the well known BDMPS-Z-W spec- 



(a) (b) 
Figure 2: The diagrammatic representation of the two first terms in (27). 
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trum, in the case of uncorrelated sources. This approximation assumes that the 
scattering centers at different times in the medium are independent^, therefore 
the medium sources can be treated as Gaussian: 



Namely, that amounts to choose the statistical weight W as follows 



W[po] = / I?bo]cxp 



71 (.T+) 



(29) 



(30) 



where n{x~^) is the medium scattering center density at the time 2"+. The 
corresponding correlator for Aq , given by (4), reads: 



where 'y{x±,yj_) can be expressed as follows: 



l[x±,y^ 



(31) 



(32) 



(2^)2 7 fci 

and is related to the dipole cross-section by the relation'' 

a{xi_ -y±) = — \'y{x^,x±_) +j{y^,yj_) + 2-/{xj_,y^)] . (33) 

In this approximation the average of two Wilson lines is a color singlet: 

(Uac{z+, z'+, x^)ul{z+, z'+,y^)) = (trt/(z+, z'+ ,x^)U\z+ , z'+ ,y ^)) 

(34) 

and the two-point function reads 



(tr[/(z+, z'+,xj_)UHz+, y^)) = exp 



d^n{^)a{x±_ - y^) 



(35) 



The second term in (27) is easily evaluated leading to (sec Appendix A for 
details): 



m - 



^ (trg.(z'+, y^; z+, x^)UHz'+,z+, z^)) 
= K^{z'+,y^ - z+ , xx_ - z±) 



^ d?(yn(0- 2^(e)a(rj.)) 



(36) 



^This approximation holds when the range of one scattering center is much smaller than 
the mean free path of the radiated gluon and of the hard parton inside the medium. 

^for more details on Wilson line averages in the gaussian approximation see ref. [47,48] 
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where r±{z'~^) = — z±^ and rj^(z+) = x±^ — z^. 

Using the foUowing property of the Green's functions, equivalent to (18), 

(37) 

and the locahty in time of the source average, the first term in (27) can be 
factorized as foUows 

j^:^ {trU{z+,z'+,z^)gliL, y^; ^^^^(L, x^;z+, z^)) 

X (Gl'^-^iL, y^;z+,u^)g'j{L, xr,z+, z^)) , 
1 

" (A^2 _ 1)2 

X {tvgl{L,y^;z+ ,u^)g^{L,xr,z+ ,z^)) . (38) 



<fu^ {trU{z+, 2'+, z^)gl{z+, u^; z'+, z'^)) 



Putting everything together, and foUowing the calculations in Appendix A, we 
recover the well-known induced radiative gluon spectrum (see for instance [30], 
where it is shown to be equivalent to Wiedemann's formulation [44]) 



d{N) 



asCR 



2SRe / dz+ / d^itj^e^'^-L "^ 



dz'+e-i ''«"(«)"("-)ai,.ai„/C^(2+, ttj.; = 0) 



^ -1 
-UJ 

-2^-.d^ylC^{L, u^; z+, y^ = 0) 
91 



(39) 



4 gluon production in proton-nucleus collisions 

In this section, we will re-derive the gluon production in proton-nucleus collisions 
in the high energy limit [45,46,49-51]. We end up by a discussion on the validity 
of the eikonal approximation at RHIC. 



4.1 Gluon production in the high energy hmit 

At high energy the nucleus is Lorentz contracted, so that we take the limit 
L —f 0, and put to = — oo. As a consequence the produced gluon is eikonal 
during the interaction time L, and the retarded gluon propagator (15) simply 
reads 

G(x, y) = h{x+ - y+)e{x~ - y-)5{x^ - y^)Uix+,y+,x^). (40) 
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Thus, the gluon field (14), when amputated of its final free propatagor, reduces 
to 



U5A\x) 

- e{x-)e{^x+)d'p{xA_) - e{x-)e{x+)d\up{x^)) 

where U = ?7(+oo, — cx); x\^. Then the Fourier transform gives 



2i 



k\ 



2{q+ + ie){q- + ie) kf 



(41) 



(42) 



where k2± ^ (1± ki± and ^pi.oton(9) is the Fourier transform of the gauge field 
of a proton alone, i.e. the Fourier transform of eq. (42) taking U = 1. The two 
terms in (42) are illustrated in fig. (3). This expression leads to the standard 
result for gluon production in proton-nucleus collisions. 




(a) (b) 

Figure 3: The two diagrams contributing to gluon production in proton- nucleus 
collisions in the limit L ^ 0. In this limit the gluon is eikonal when passing 
through the nucleus, and the gluon emission inside the nucleus is neglected. 



4.2 Validity of the eikonal approximation 

The center of mass energy per nucleon at RHIC is 200 GeV. It corresponds to 
a Lorentz contraction factor of about 7 = v^/2mp ~ 100. For a nuclear radius 
Ra — 6.5 fm, we end up with the estimate L ~ 0.5 GeV^^. For the eikonal 
approximation to be valid, the gluon production time has to be much bigger 
than the maximum interaction time: 

tprod - > L. (43) 
91 
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At mid-rapidity, at RHIC, uj ^ q± ^ Qs ^ IGeV. We get tprod of order of L. 
Therefore, one probably has to go beyond the eikonal approximation for the 
gluon propagator to describe mi-rapidity data at RHIC. 

The corrections to the eikonal approximation are taken into account in formula 
(27) for gluon radiation. Even if (27) has been derived in the framework of final 
state interactions in nucleus- nucleus, it is applicable for gluon production in 
proton-nucleus collisions (assuming the medium to be a nucleus and the projec- 
tile to be a proton, and taking the projectile production time to = — oo). This 
is a possible phenomenological application for mid-rapidity hadron production 
at RHIC. 

At forward rapidity, uj ~ cj±e^, therefore, uj is enhanced by a large factor (ex- 
ponential of the rapidity) leading to a much larger gluon production time com- 
patible with the eikonal approximation. 



5 The relation between q and Q. 



The transport coefRcient q which characterizes the density of the scattering 
centers (gluons in the medium) in the medium produced in nucleus-nucleus 
collisions is defined to logarithmic accuracy by ^q{£,)'''\ — '^(0'^(''-l) [53]. Using 
eq. (35), we can write the following relation: 

^J—{trU\x+,y+;z^)U{x+,y+;z'^))=c^p J^^ dm)i^±-^'i^) 

(44) 

Whereas q is local in time, the saturation scale Qs is a global quantity defined 
as [48]: 



-(trC/t(i,0;2i)[/(i,0;2l)) =cxp 



1 



(45) 



where L = 2Ra is the nuclear diameter. From this formal analogy it becomes 
clear that the saturation scale appears as an initial condition for the transport 
coefficient in nucleus-nucleus collisions as suggested in [54] 



9(e = 0) ~ qcoM - QI/2Ra- 



(46) 



Obviously the transport coefficient increases due to the strong interactions of the 
freed gluons, from cold to hot matter. This mechanism is not discussed in this 
work, however, the medium dynamics could be encoded in the time evolution 
of the medium field, which could lead eventually to thcrmalization. 



6 Summary 

In this paper, we calculate the gluon spectrum of a high energy and point-like 
projectile passing through a dense medium, presenting a compact and simple 
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derivation based on the solution of the Yang-Mills equations in the light cone 
gauge of the fast projectile. We reproduce in a straightforward way the gluon 
spectrum in proton-nucleus collisions at high energy and the induced gluon spec- 
trum in the final state of nucleus- nucleus collisions. We also discuss the validity 
of the eikonal approximation in gluon production in proton-nucleus collisions at 
RHIC, and argue that one could have non-negligible non-cikonal contributions 
at mid-rapidity gluon production. As an application of our formula for gluon 
production beyond the eikonal approximation (27) to phenomenology, it could 
be interesting to see whether this complete solution leads to an appreciably dif- 
ferent answer for hadron production spectra and for the Cronin peak at RHIC. 
Finally, we recall the relation between the saturation scale in initial state inter- 
actions and the transport coefficient in final state interactions, as the transition 
from cold to hot matter in nucleus- nucleus collisions. 
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A Averages over medium sources 



In this appendix we shall derive analytic expressions for some useful averages 
over the medium sources in the Gaussian approximation. 

A = J^^^ {tYg^{x+,x^; y+, y ^)U\x+ , y+ ; z ±)) 



'Dr±{^) exp 



y 



X {trU{x+, y+; r^(0)C/t(x+, y+; z^)) , (A-1) 

where r±{x^) = x±_ and t"j_(j/+) = yj^. In the gaussian approximation the 
two-point function can be easily evaluated and gives 

1 



j;^^ {trUix+,y+;r^{0)UHx+, y+; z^)) 



exp 



(A-2) 



Now we define 



ICUx+,xr,y+,y^) = Jvr^iOcxp 1^ d^(|^fl(0-in(e)a(rx)) 



(A-3) 



13 



Making v±^ = r± — zj_, we get 

j^2^_ ^ {trgu,{x+,x±;y+,yj_)U''{x+ ,y+;z±)) = ICuj{x+ ,x±-z±;y+ ,y^-z±). 

(A-4) 

Now let evaluate the following useful average 
1 



B 



- {trgl{x+,xr,y+, y^)g^{x+,zr,y+, z'^)) 



Vr^iO / I?Mi(e)exp 



lUJ 

T 



: {trUHx+,y+; r^)U{x+,y+; u^)) 



(A-5) 



with the following boundary conditions: r±{x^) = r±{y^) = y^, it^(x+) = 
z± and u±{y^) = z'j_. Recalling that the two-point function depends only on 
rj_ — we perform the following change of variables 



oi.±_ = It J 
(3± = uj 



- rj 

- r I 



yielding 



B 



Now we can perform the /3 integration yielding to a (5-function constraining the 
ct variable to a straight line as following 



1 



[(e - y+)a±{x+) + ix+ - C)ax(jy+)] • (A-7) 



It is straightforward to show, for instance by discretizing the path integral, that 



B 



2'k{x — y)+ 
iuj 



X exp 



— —ilx - y)i -(z- z')\) - i 



We shall calculate the Fourier transform appearing in (27) 



(A-8) 



(A-9) 
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Let 



vj_ = (x - y - z + z')±, 
w± ^ {x - y + z - z')±, 



(A-10) 



Wc act 



OJ 



2t:{x — y) + 
X exp 



1 r 



where now 



^±(0= P ^2 ^± + {y-z') 



(A-ll) 



(A-12) 



(x - y)+ 

o'{l±) is independent of w±, thus the integral over w± yields a (5(t'_L) and we 
finally get 



C = cxp 



, 1 

-ikj_.{y - z')j_ - - J ^ d^n{^)a{{y - z')j_) 



(A-13) 
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